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The kinetics of an initially undercooled solid-liquid melt is studied by means of a generalized 
Phase Field model, which describes the dynamics of an ordering non-conserved field 4> (^-g- solid- 
' liquid order parameter) coupled to a conserved field (e.g. thermal field). After obtaining the rules 

^\ , governing the evolution process, by means of analytical arguments, we present a discussion of the 

0^ ' asymptotic time-dependent solutions. The full solutions of the exact self-consistent equations for 

the model are also obtained and compared with computer simulation results. In addition, in order 
f— I , to check the validity of the present model we confronted its predictions against those of the standard 

' Phase field mode and found reasonable agreement. Interestingly, we find that the system relaxes 

towards a mixed phase, depending on the average value of the conserved field, i.e. on the initial 
condition. Such a phase is characterized by large fluctuations of the field. 



PACS numbers: 64.60C, 64.60M, 64.60A 



^ ■ In the last few years considerable effort has been devoted to the study of systems far from equilibrium Well 
' known examples are provided by phase separating systems, initially prepared in a state of equilibrium, and rendered 
\l , unstable by modifying a control parameter such as temperature, pressure, or magnetic field. To restore stability they 
evolve towards a different equilibrium state determined by the final value of the controlling fields. Such evolution can 
, be very slow and is often characterized by non uniform, complex structures both in space and time. 
— . ■ Two simple dynamical models, often called Model A [|| and model B have been introduced in the literature in 
^\ \ order to understand kinetic ordering phenomena (see also [^D). The model A describes the growth process, when 
"Y^ ' the order parameter is non conserved, whereas the model B is appropriate if the order parameter is conserved. In 
the first case, the late stage growth is driven by the tendency of the system to minimize the energy cost due to the 
^ . presence of interfaces between regions separating different phases. Thus, as the curvature decreases the process slows 
down and the domain size, L{t), grows in time according to the law Ht) ~ t^/"^. In the conserved case, instead, the 
approach to equilibrium is limited by the diffusion of the the aggregating material, as larger domains can grow only 
at the expenses of smaller ones. The average size L[t) increases proportionally to t^/^, where the dynamical exponent 
^ ' z is 3 for scalar order parameters and 4 for vector order parameters. 

A further model, known as Phase field model (PFM), is somehow intermediate between A and B and consists of two 
fields coupled bilinearly: one field represents a non conserved ordering parameter, with type A dynamics, whereas the 
second is a temperature shift field subject to a diffusion equation supplemented by a source term. The model can be 
?H ' cast in the form of coupled partial differential equations for a non conserved order parameter interacting with a time 
dependent conserved field. Its dynamics is very rich, since it displays features characterizing both the pure A and the 
pure B models as it is revealed from the analysis of the structure functions at different times. In other words, after a 
rapid initial evolution one observes an intermediate stage in which the growth is curvature driven and an asymptotic 
regime during which diffusion limited behavior is seen. 

The PFM, introduced and physically motivated by Langer provides a theoretical framework for many natural 
processes. It is designed to treat situations where the relaxation dynamics of the order parameter associated with the 
presence of a liquid or a solid is coupled to the diffusion of heat released during the change of state. 

An example is the growth of a solid nucleus from its undercooled melt, a phenomenon encountered in rapidly 
solidifying materials, such as metals, where the growth is limited by the rate of transport of the heat of fusion away 
from the solid- liquid boundary [Q. As the heat released by the solid accumulates at the interface, it slows down 
the growth, because diffusion must act over a thicker and thicker region. This mechanism has also implications in 
the morphology of the growing phases and is responsible for the instability of a planar solid-liquid interface with 
respect to a perturbation of its shape; one realizes, immediately, that a protrusion of the solid phase into the liquid 
advances faster than its neighboring regions, because it explores a region where the undercooling is greater, so its 
growth becomes faster. The solid-liquid surface tension eventually provides the necessary balance and prevents the 
interface to be eroded by fluctuations of very short wavelength [||,^. 
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Another very closely related problem, is the the growth of the solid phase in multicomponent solutions, where one 
of the components is to be diffused away from the interface in order to form a stable crystal 0. In the present 
paper, we shall confine the discussion to the thermal case for the sake of clarity, and investigate the kinetic ordering 
of a spherical version of the Phase field model. This study, extends our previous investigations llo|-[l^ to include 
a non vanishing order parameter. In our opinion it can be useful, because it provides one the few models, whose 
static properties can be obtained exactly in arbitrary dimensionality and whose relaxation behavior can be analysed 
in great detail by means of analytical and numerical methods. We note that in the field of ordering kinetics there 
exist only a few models for which the relaxation can be studied without performing heavy numerical calculations. 
In particular the late stage behavior of processes with conserved dynamics is hardly observable numerically, due to 
computer limitations ||l^,|l3|. Besides the examples cited above, this study may be of some help to treat analytically 
some models introduced recently with the aim of describing irreversible aggregation phenomena [ p^ . 

In the present paper we generalize the model, introduced previosly by two of us to the case of off-critical 

quenches, e.g. to initial conditions corresponding to non-vanishing values of the fields. The structure of the paper is 
the following. In order to make the paper self-contained we have included two sections (I and II), where we recall some 
basic notions, which lead to the thermodynamic derivation of the Phase field model ||6|,|l5|-|l7| . and the construction 
of the Lyapunov functional, from which the coupled equations of motion of the PFM model can be derived . In 
section III, we state the Spherical Phase field model and write explicitly the closed set of equations, which we discuss 
qualitatively in section IV . In section V comparisons with numerical simulations in d = 2 at zero temperature are 
illustrated. The predictions of the spherical model are confronted with those of a more realistic scalar order parameter 
Phase field model and the similarities and differences are stressed in the conclusions. 



I. THERMODYNAMIC PRELIMINARIES 



Let us consider a material which under suitable conditions of pressure and temperature can exist in two distinct 
thermodynamic phases, a liquid and a solid. If the pressure is held constant at the value corresponding to solid-liquid 
equilibrium and the temperature, T, is varied one can favor the solid phase for T < Tm, or the liquid phase for 
T > Tm, where Tm is the melting temperature at which the equilibrium first order transition takes place. One 
usually calls undercooled melt a material brought below its melting temperature, but still in the liquid phase. 

We shall consider the situation T < Tm, which is experimentally and technologically more interesting. Below Tm 
the value of the thermodynamic Gibbs potential of the solid phase is lower than that of the liquid phase, which is only 
metastable. A convenient way of studying the solidification process is to adopt a phenomenological Ginzburg-Landau 
description by introducing a suitable crystalline order parameter <j) which assumes the conventional value 4>i in the 
liquid phase and (ps in the solid phase. One then employs a field theoretic free energy functional of the form: 



F[(j)\ 



^(V</.)2 + /(0) (1) 



where /(0) is a function of the order parameter with the property of having two minima of equal depth a,i (j) — (pi 
and 4> = (j)s- 

The constant e has dimensions of energy/volume and is for the moment arbitrary. The gradient term represents 
the energy cost necessary to create an inhomogeneity in the system, the quantity ^ has the dimension of length and 
is associated with the scale over which a inhomogeneity in the system vanishes. Upon minimizing F[<f>\ with respect 
to (j) and selecting the non uniform localized solution of the variational Euler-Lagrange equation corresponding to the 
lowest value of the Gibbs free energy one obtains the surface tension a of the model which is proportional to the 
correlation length 

a ^ (2) 

and thus to the interface thickness. The numerical coefficient is of order 1 and will be ignored, because it does not 
influence our discussion. 

In order to include undercooling or overheating effects, i.e., a temperature different from Tm, we introduce a 
dimensionless field 

u{x)^^j-[T{x)-TM] (3) 

proportional to the local temperature shift (T(x) — Tm)- The constants Cp and L are respectively the specific heat at 
constant pressure and the latent heat of fusion per unit volume. The local field u acts as an external field, conjugate 
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to the crystalline order parameter (jj, favoring the solid phase for u < and the liquid phase for u > 0. For u = the 
two phases coexist. 

As usually done for first order phase transition, metastability is taken into account by eliminating (p in favour of u 
via a Legendre transform. One then introduces the Gibbs potential 

G[u] = F[(j)] -XeJ d'^xu{x) 0(a;) (4) 

where (j) = ^[w] is obtained form 

and A is a nondimensional parameter. A convenient way to relate A to the known thermodynamic parameters is to 
consider the entropy difference between the pure uniform solid (0 = 0s) and liquid (0 = (pi) phases at the melting 
temperature. This is related to the latent heat by the relation 

Si-S,^^ (6) 

where Si and Sg are the entropies of the liquid phase and of the solid phase, respectively, and V the volume of the 
system. By using the thermodynamic relation 

dT- ^' 



equation and the expression (Q) for uniform fields, we get 



LV _ dG 



dG 
df 



from which we obtain: 



= XeV{<Pi-<Ps)^ (8) 
A= (9) 

where A0 = — 0s. 

Next consider a solid spherical drop of radius ^ immersed in an undercooled melt (u < 0). The Gibbs potential 
G with the drop is - see eq. ^ ~: 

G = Go + Xeu—R^A(l) + 4TrR^a (10) 
o 

where the first term is the Gibbs potential G without the droplet, the second term is gain in replacing the liquid with 
solid in the droplet, finally the third term is the cost in creating a surface separating the liquid and solid phases. In 
equilibrium no energy is needed to create the droplet so G is stationary with respect to variation of R. By imposing 
SG = we readily obtain the critical nucleation radius Rn 

Rn^^ (11) 



where 



is a capillarity length, which using the expression of A, eq.(p|), can be written as: 

2aCpTM 

do = — j2 — ■ (13) 
Finally from eqs. and (12) it follows that we can write the dimensionless parameter A as the ratio two length: 

A = -f (14) 
do 

where we have defined = (2/A0) a [cfr. eq (^)]. We see that A is small provided the interfacial thickness is much 
shorter than the capillary length. 
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II. THE PHASE FIELD MODEL 



In this section we shall introduce relaxational dynamics into the model. A large body of work in the area of dynamic 
phase transitions has been focused on the time-dependent Ginzburg-Landau (TDGL) model, because of its capability 
of describing a variety of problems. In equilibrium the field (t>{x) minimizes the Gibbs potential G. Thus we assume 
that the approach to equilibrium is described by following equation: 



d(j){x, t) 



■G[0,u] 



dt ''S(j){x,t) 



(15) 



it 



where the last equality is obtained using eg. (^. If the field u varies on time-scales much longer than those of 
can be considered "quenched" and eq. ( |l5| ) would be the standard non conserved TDGL equation, or model A. 

In the Phase Field model, and in the absence of external sources, u{x,t) is assumed to evolve on time-scales of 
the same order of magnitude as those of (f> towards an homogeneous configurations. The time-evolution of u is now 
coupled to that of </> and cannot be neglected anymore. In fact, when a piece of material solidifies it expels some heat 
and the surrounding liquid melt warms up, causing the average temperature to increase. In turn when a region of 
solid melts it adsorbs some heat is adsorbed and the liquid becomes colder. 

As a consequence eq. (|l5|) has to be supplemented with an equation for u. The thermal field u{x, t) is subject to the 
Fourier equation of diffusion of heat plus an additional source term which represents the latent heat of solidification, 
accompanying the appearance of the solid phase. 

The energy balance requires that the latent heat released at the transition equates the temperature change of the 
melt multiplied the specific heat, i.e: 



du{x, t) 
dt 



DV^u{x,t) 



1 d(f>ix,t) 
dt 



(16) 



where D is the thermal diffusivity and the last term in the right-hand side is the amount of material which crystallises 
per unit time and thus proportional to the heat released during the first order transition |^. The coefficient Atf> 
guarantees the correct energy balance. Notice that the last term represents a source of heat when d(j>{x,t)/dt is 
negative, i.e. when the system solidifies, or a sink when it melts, positive d(t>{x,t)/dt. In other words, since we 
are considering a closed system the total amount of solid produced is proportional to the change of the average 
temperature of the system. 

The two dynamical equations ( p^ and ( p^ can be obtained from a unique Lyapounov functional J-, which plays 
the role of the time dependent Ginzburg-Landau potential in the present problem. In order to establish the form of 
J- we perform the transformation 



A(j) 



and eliminates u in favor of the new field U. One can then write eqs. (|l|) and (|l|) as 



d(l>{x,t) 
dU{x,t) 



5T 



dt 



6(t){x, t) 
ST 



dU{x,t) 



with the Lyapounov functional | fi3| , p2[ 

U] = / d'^x 



A0 



(17) 

(18) 
(19) 

(20) 



Note that the dynamics of U is conserved. When the temperature field vanishes, i.e., U = 4>/A(f>, the functional J" 
has two equivalent minima, corresponding to two spatially uniform solutions: the uniform solid and liquid phases. In 
general eqs. (p^-(|20|) generate a complex dynamical behavior which has been the object of some studies. 

In the long time limit we may expect that while u becomes homogeneous, the crystalline field (j) roughly assumes 
only the two values and (j)s- If this is the case, from the knowledge of U and u we can compute the fraction of 
volume occupied by the two phases. Indeed we can write (p — Xs4>s + where Xs and xi are the fraction of volume 
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occupied by the solid and liquid phase, respectively, and the overbar denotes spatial average. From eq. 
condition Xg + xi = 1, we get 



7f) and the 



Xl 



7^ 



U 



(21) 



where u ~ u is the asymptotic value, while U is the initial value being its dynamic conserved. From this it follows that 
if the asymptotic value of u is zero, i.e., the system relaxes towards a two phase coexistence, the fraction of volume 
occupied by each phase is determined only by the initial value of U . We also notice that when the system starts with 
an undercooling u = — 1, and order parameter m{t = 0) = 1, i.e. C/, the latent heat produced is just enough to heat 
the melt at the final equilibrium temperature, u — Q. In such a case the final volume fraction of the solid is simply 
Xs = 1/2 — ?7, and attains its maximum for U — —1/2. 

The functional T decreases with time as it can be shown using the equations ( |l^ ) and (^9|): 



'dt 



5T d(j){x,t) 5T dU{x,t) 



ix,t) 



- / d'^a 



dt 



SU{x,t) dt 



{x,t) 



D 



VST 
SU{x,t) 



< 0. 



(22) 



So far we have discussed purely deterministic evolution of the order parameter and of the thermal field. Noises can 
be added to both equations to represent the effect of short wavelength fluctuations; in this case eq. ( p2|) does not 
hold. 



III. SPHERICAL PHASE FIELD MODEL 



The choice of the local function f{(f>) is somehow arbitrary, as long as the general property 

lim f{4>) = +00, and two equal minima for (p = (ps and (p — (pi (23) 

are satisfied. A largely used form for f{(p) is 

m) - 9 - ^) (24) 

which is even and has two equal minima aX <p = (ps ^ —1 and cp ~ cpi — 1^ so that A(/) = 2. The parameter g gives the 
strength of the local constraint. In the limit of large positive value of g the field cp can take only the values (ps and cpi 
(Ising-like variables) . 

The phase field model described by equations (|l8|)-(^^ contains all the relevant ingredients necessary to describe 
the phase separation occurring in solid forming melts. However, due to the local nonlinear terms contained in the 
function /(0) the solution of the dynamical equations (p^-(pO|) are far too difficult and are known only for some 
special situations. In the general case the known results follow from dimensional arguments. 

To overcome this difficulty, an alternative strategy is to modify the model into a simpler one, yet maintaining the 
general properties. This can be achieved by replacing the local quartic term in (p4): 



d<'xd^{xf^^ 



d X (pix 



(25) 



This kind of constraint is much softer than (^4|) since it does not act on each site, b ut g lobally over the whole volume. 
In the following we shall denote it as globally constrained model or spherical model jTsp^ in contrast with the model 
(p^)-(pO|) where the constraint is local pl| , p4 |. 

The price one pays for this change is the loss of sharp interfaces between two coexisting phases. As Abraham and 
Robert j2^ showed longtime ago the spherical model in zero external field displays two ordered phases below the critical 
temperature, but no phase separation. Equivalently one can say that a planar interface between two coexisting phases 
is unstable, due to the presence of long-wavelength excitations analogous to spin-waves, an instability much stronger 
than the one due to the presence of capillary waves in the scalar order parameter case p4,M. As a consequence. 
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while this choice is very convenient for analytic calculations, it changes the structures of the non uniform solutions in 
the static limit. Nevertheless, in spite of this fact, the model has a rich phenomenology as we shall see below and the 
approach to equilibrium remains highly non trivial. 

From equations (po|), (|2^ ) and (^5|), the potential T for the spherical model reads: 



which substituted into: 



d(j){x, t) 
d1~ 



X 



g] r + AC/^ - XU(f> 



U]+i^{x,t) 



{x,t) 
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(26) 



(27) 



dU(x,t) 
dt 



(28) 



determines the time evolution of the fields and U . We added to the evolution equations a noise term to simulate 
the effect of short wavelength fluctuations. The two fields 77 and f are independent Gaussian fields with zero mean 
and two-point correlations: 



(7?(a;, t) rj{x', t')) =2TfT^ 6{x - x') S{t - t') 
{C{x, t) ^{x', t')) ^-2TfD V^S{x - x') 5{t - t') 
{r^ix,t)ax',t'))^0 



(29) 
(30) 
(31) 



where Tf is the temperature of the final equilibrium state whereas D and are the kinetic coefficients appearing 
into eqs. ( p7| ) and (|2^). 

It is useful to separate out the spatially uniform component of fields (f> and U. Thus introducing the Fourier 
component of the fields we have: 



E 

k 



(j){k)e 



ik-x 



(X e 



— ik-: 



(32) 



and 



0(fc,i) 
U{k,t) 



Q{t)5k,o - 



-S(j){k,t) 
SU{k,t) 



(33) 
(34) 



where both 5(j) and SU are zero for k = \k\ = 0, ■m{t) ~ (f>{t), Q{t) = U{t). 

To study the behaviour at finite temperature Tf it is also useful to introduce the equations of motion for the three 
equal-time real space connected correlation functions C^^{r, t) = {(f>{R + r, t)(f){R, t))c, C^u{f, t) = {(f){R+r, t)U{R, t))c 
and C(ju{r, t) — {U{R + r, t)U{R, t))c, whose Fourier transforms are the structure functions. The average <> is over 
the external noises r/ and ^ and initial conditions. 

Due to the special form of the non-linear term in the equation of motion the set of evolution equations for the 
averages m{t), U{t) and the correlation functions is closed. Indeed, in the Fourier space these read: 



d(j){k,t) 

dt 
dU{k,t) 

di 



F^{k)+r]{k,t) 
Fu{k)+ak,t) 



(35) 
(36) 



where i^,/,,;/ are the Fourier transforms of the first term on the right-hand sides of eqs. (| 
have 



and (28). From iM) we 



F^{k) = M^^k, t) cp{k, t) + M^u{k, t) U{k, t) 
Fuik) ^ Mu4k,t)(j){k,t) + Muuik,t)Uik,t). 



(37) 
(38) 
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where the matrix elements are given by 



M^^(k, t) = -T^iek"^ + r- + gm^{t) + gS{t)], 

Mu4>{k,t)= D\k'\ 
Muuik,t) = -2DXk^. 



(39) 



where r — {—g + A/2) and, in the hmit ^ oo the quantity S(t) is the integrated (/)-structure function 



= J2Wk,t)^{-k,t)), (40) 



and m(i) — (</>). 

As a consequence we have 



dm{t) 

dt 
dQ{t) 



M^^O, t)m{t) + M^t/(0, t)Q{t) (41) 
(42) 



dt 

expressing the conserved nature of the field J7, and 

1 d 

--g-^C^4k,t) = M^4k,t)C^^{k,t) 

+ M^uik,t)C^uik,t)+T^Tf (43) 

— C^uik,t) = Mu4k,t)C^4k,t) 

+ [Muuik,t) + M^4,{k,t)]C4,u{k,t) 

+ M^u{k,t)Cuu{k,t) (44) 

1 d 

2 ^ Cuu{k, t) = Mu^{k, t)C^u{k, t) 

+ Muu{k,t)Cuu{k,t)+DT}e. (45) 

We note that a closure at the same level would have been obtained in the framework of a Hartree approximation 
for the model with local constraint described by eqs. (^-(^. However, within the present model the eqs. (pll)-(|45|) 
are exact and not the result of a approximate decoupling of the correlations. 



IV. LONG TIME BEHAVIOUR 



In this section we shall discuss the behaviour of the spherical Phase Field model for long times, i.e., t ^ oo. The 
results will be compared with those of direct numerical simulation in the next section. 

We assume that at the initial time we have an undercooled liquid with some supercritical solid seeds. This means 
that M < while m = (j) lies in the interval (0, 1). For a generic initial configuration, the undercooled liquid is not in 
equilibrium. Thus at the initial stage the relaxation of the field (f) is only slightly modified by the dynamics of the 
slower field U , which can be considered "almost quenched". During this stage the size of the solid seeds grow with 
time, while the maximum of the structure function is located at fc = and grows with time. If this regime is long 
enough, one can recognize a typical non conserved order parameter dynamics (NCOP) domain growth proportional 

to 

This kind of behaviour persists until the typical size of the domains reaches that associated to the conserved U field. 
At this time the dynamics of the two fields becomes strongly correlated, and the conserved order parameter dynamics 
eventually dominates. As a consequence the <j) field slows down, since the coupling with the conserved U field acts as 
an additional constraint, while m{t) becomes nearly constant, and equal to the asymptotic value, i.e., m ^ <p^. 
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The cross-over time tc can be readily estimated from dimensional arguments. An inspection of equations (|35|)-(|3q) 
reveals that under a suitable transformation of parameters, in which Xt —>■ t, the dimensionless parameter A can be 
traced out. 

This means that the crossover time, tc, to the conserved dynamics is also of order ~ 1/A. It can be shown that the 



Langevin equations (35)-(p6|) obey detailed balance and that the stationary probability density is 



Tf 

Using this we easily get the equilibrium (pcj) correlation function, which reads 

Tf 

^^k^ + r + gS + gm^ - X/2' 



Pst[<P,U]^exp{~^T[cl>,U]] . (46) 



(47) 



The appearance of an ordered phase, with m ^ for temperatures Tf, below the critical temperature Tc, is revealed 
by the divergence of the fc = mode. This implies that the equilibrium ( i.e. when t oo) value of m must satisfy 
the following equation: 

r + gS + gm^ - A/2 = 0. (48) 
On the other hand, from the equation of motion ( ^l| ) m must be solution of: 

rm + girC' + gSm - AQ = 0. (49) 
The simultaneous solution of eqs. (^8|)-(^9|) requires: 

m = 2Q (50) 

which via eq. ( ^7|) together with — 2 implies u = 0. This means that the system relaxes towards a nontrivial 
phase coexistence state, with non vanishing order parameter and diverging small k fluctuations. 

One sees that the condition ( |50| ) can be satisfied only for — 1/2 < Q < 1/2. If Q lies outside this interval the system 
does not relax to a mixed phase, but instead settles in a spatially uniform state without zero modes. Indeed, in this 
case S{t) vanishes and m relaxes for long times to the value given by, cfr. (p]), 

rm + gm^ — XQ = 

~gm ( 1 - m^) - Au = 0. (51) 

We note that eq. (^ij) is equivalent to say that a a spatially uniform field (j>{x) is a stable minimum of the the potential 
T. For values large enough of equation ( ^ ) has a single solution, positive for u > and negative for u < 0. These 
correspond to a liquid phase above the melting temperature, i.e., positive m and u, and a solid phase below the melting 
temperature, negative m and u. As decreases two additional solutions eventually appear. One of these is unstable 
while the other represents a metastable spatially uniform state: i.e. solid above the melting temperature {u > 0) and 
liquid below the melting temperature (m < 0) . These solutions, however, are unstable against fluctuations, indeed the 
presence of S{t), which vanishes only for t oo prevents the dynamics to reach these metastable states. Therefore 
for any value of u the physical solution of eq. ( pl| ) is the most negative one for It < and the most positive one for 
u > 0. For small u we have m = (1 — 2{X/g)u) signu. 

A more detailed analysis of the approach to equilibrium can be done employing a quasilinearization procedure, i.e. 
we assume that the quantity R{t) — r + gm{t)^ + gS{t) can be treated as a constant, along different pieces of the 
trajectory |]l0| . One can verify at posteriori that the assumption is valid and leads to useful predictions. Since the 
behaviour at = is representative of the entire dynamics in the ordered phase when Tf < we also set Tf = 0, 
without loosing relevant information. 

Assuming the quantity R{t) to be nearly constant eqs. (|3^) and ( ^6| ) become a linear system whose solution has 
the form: 

0(fc,t) =c+(fc)e-+W* + c^(A;)e--W* 

U{k,t) = c+(fc)e'^+W* +c[}(fc)e'^-('=)* ^ ' 



where w+(A) and w_(fc) arc the eigenvalues of the M matrix 



T^i^'^k^ + R)- 2DAfc2 ± /[r^(^2fc2 + R) + 2DAfc2]2 + AT^DX^k^ 



(53) 
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For time t ^ 1 the dynamical behavior of the solution is determined by the larger eigenvalue uj+ (k) For large values 
of the eigenvalue uj+(k) decreases as —k^, thus to discuss the behavior of the solutions after the initial transient 
a small k expansion of uj+{k) is sufficient. The form of this expansion depends on the sign of R{t). When R{t) is 
negative the appropriate expansion of i-Uj^{k) is: 



^{k)=T^\R\- 



(54) 



Notice that in thie regime where eq. is valid there is a competition between the curvature term F^^^, which 

represents the driving force of the dynamics of the pure Model A, and the term DA^/ji?] due to the coupling to 
the heat diffusion. For R{t) > 0, the representation eq. (|5^ ) breaks down and one one must instead consider the 
expansion: 

u+{k) = 2DX (^-^ -l^k^- ak'^ (55) 

where C4 is a positive coefficient. 

Since the value of R{t) as well as its sign change along the trajectory, one must employ either eq. (Q) or eq. ( [55| ) 
depending on the stage of the growth process. One the other hand, this kind of analysis is not applicable in the 
crossover region where m{t) and S{t) vary too rapidly and R{t) ^ 0, but this fact does not invalidate our findings. 

All the relevant behaviors can be classified according to the value of the conserved field Q\ to this purpose one must 
treat separately the cases Q > and Q < 0. 

Let us consider first the case Q < and assume ra{t = 0) ^ 1, i.e. the system initially is formed by an undercooled 
liquid. In the initial stage the fields (j) and U are nearly uniform and characterized by small fluctuations. In the 
evolution equation for m , U plays the role of a constant field of value Q, hence 

— = -F^ [Rit)m{t) - XQ] (56) 

describes the relaxation of m in a static field Q. During such a stage the system relaxes towards the nearest fixed 
point which makes the right-hand side of ( ^6|) to vanish. Equating to zero dm/dt we get the relation 

R{t) = (57) 
m 

which is negative for m > and Q < 0, so that the relevant expansion is (|54|). Thus the liquid phase is unstable as it 
can be seen from eq. (^4[); in fact, as long as m{t) remains positive the system develops strong fluctuations about the 
uniform mode, fc = 0. Such a regime lasts for a time of order ~ 1/A, after which the fluctuations eventually drive 
the system towards negative values of the order parameter as to reduce the free energy cost. 

In the successive stage the system becomes prevalently solid, and this is signalled by the change of sign of m. Also 
R changes sign so the expansions (jsj) becomes invalid. However, for times t ^ 1/A, long after the transition, the 
evolution of m slows down again, i.e. dm/dt ~ and i? > 0, so that ( p5|) is appropriate. Using the result eq.(|57|) the 
relevant form of w_|_(fc) is: 

Lo+{k) = Dx(^^-l^k^-Cik\ (58) 

Such expression shows that whenever m/2Q < 1 all modes of flnite wavevector are damped and one cannot observe 
growing modes: the system relaxes towards a spatially uniform state. 

On the other hand if m/2Q > 1 fluctuations are large up to a finite wavevector, because the fastest growing mode is 
located at a finite k. This case corresponds to a phase spatially non uniform with large fluctuations. Asymptotically 
m = 2Q and the peak position km, moves towards vanishing values of k. This scenario is typical of the conserved 
order parameter (COP) dynamics. Indeed, as done in Ref. |ic| , [ll|] , it can be shown that in this regime the dynamics 
exhibits multiscaling ||lO| . 

In Fig. |l| we report the behavior of dm/dt versus the order parameter m. For small negative values of XQ, dm/dt 
nearly touches the horizontal axis, and the evolution of m becomes very slow. 

Let us now turn to the Q > case. In this case from ( p7| ) follows that R{t) > 0, so that the relevant expansion is 
(p8|). In analogy with the Q < case, from (|58| ) we conclude that if Q > 1/2, the phase to 1 is stable and the final 
state is pure, i.e. is characterized by small spatial fluctuations. 

For smaller value of Q, i.e., < Q < 1/2 the coefficient of the term of order k^ of the eigenvalue uj+ is positive and 
the peak of C^^{k,t) is located at a finite wavevector. This causes an instability of the initial pure phase to 1 and 
the appearance of a mixed phase, characterized by a lower, but still positive, value of m{t). 
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V. NUMERICAL RESULTS 



In what foUows we shaU consider the zero temperature dynamics, because it is known to be representative of the 
subcritical behavior. We have compared the numerical results with the predictions of the previous section and found 
good agreement. 

We have solved numerically the equations of motion (^5|)-(^6|) in two dimensions by using a simple Euler second 
order algorithm and a discretization of the integrals on a N x N bidimensional lattice. We used periodic boundary 
conditions. The parameters employed are 

A = 4A/e2 = i/a (59) 

D = e2/8A 

with e = 0.005, A = 0.1 and a = 10. The time step used is At = 2 x 10^^ and a lattice space Sx = 0.01. Different 
values of N were used, here we report the results for N = 256. 

The system was initially prepared in an non uniform initial state formed by 70% of undercooled liquid and the 
remaining 30% of seeds of solid randomly distributed. This ensures that at the initial time the order parameter 
m{t = 0) is positive and the correlations are small. As long as these two conditions are met, the results are not too 
sensitive to the initial solid/liquid fraction. The thermal field u was taken uniform, Uij(t = 0) = u(t = 0), with both 
positive and negative values. According to eq. (U^), the conserved field U was chosen to be equal to 



U^,{t = 0) = u(i = 0) + - 0y(< = 0). (60) 



1 
2 

The model exhibits two different long-time regimes as far as the temperature field is concerned. The snapshots of 
the field indeed reveal that changing from positive values of the conserved field Q to negative values the long-time 
morphology changes. For Q > one observes solid drops immersed in a liquid and u < 0, whereas for Q < liquid 
drops become trapped in a solid matrix and u > 0. This is in agreement with the general results of thermodynamics. 
For quenches with positive Q the minority phase is the solid, thus we expect the solid drops to have positive curvature. 
Indeed, relating the curvature /C, to the temperature field we obtain from the Gibbs-Thompson condition u = —dglC 
a negative value of u for Q > 0. On the contrary, upon crossing the Q = line, the solid becomes the majority 
phase and u change sign, since the curvature relative to the solid is negative. 

In Fig. 1^ we show the behaviour of the order parameter m as function of time for the four distinct regimes. In all 
cases at the initial time we have m = 0.54. We note that in the case \Q\ > 1/2 the system evolves towards a spatially 
uniform state, which is a liquid above the melting temperature for Q > 1/2 and positive initial u(t = 0), case (a) in 
figure, or a solid below the melting temperature for Q < 1/2 and negative initial u{t = 0), case (d) in figure. This 
state minimizes the potential 

In the case \Q\ < 1/2 the system evolves towards a phase equilibrium state at the melting temperature: we have 
for t^oou = u = and m — 2Q, in agreement with the analytical results of the previous section. 

When the systems starts from an undercooled state, i.e. u < at the initial time - cases (b), (c) and (d) in Fig. || - 
despite the fact the liquid state is unfavorable, m increases for short times. The system then evolves towards a liquid 
state, m tends to saturate to a fixed value ~ 1 and the fiuctuations are small. This state is however unstable with 
respect to fiuctuations, and indeed after a time tc we observe a transition towards the asymptotically stable state, 
which depends on the value of Q < 1/2. 

The time is a decreasing function of A, as can be seen in Fig. |3[ where we report m as a function of time for 
— 1/2<(5<0 and variuous values of A. In Fig. ^ we report tc versus A for the curves of Fig. ||. The Hne is the 
theoretical prediction ~ 1/A. In this figure we defined tc as the time such that mint Tn{t) = 2Q. Other definitions 
are possible, for example when m(t) = 0, or any other fixed value. All these definitions leads to the same scaling with 
A. 

The scenario here depicted remains valid also for the local constraint case. This has been checked by numerically 
integrating the appropriate equations. The main difference between the local and the global case shows up at mor- 
phological level, as can be seen in the snapshots of Figs ^ and ^ As one can see while the phase field model with local 
constraint has sharp domain wall the model with the global constraint presents smoother interfaces. We stress that 
despite this difference the circularly averaged correlation functions and the structure factors for two cases are quite 
similar. In Figs. |^ (global) and ^ (local) we report the circularly averaged cj)4> correlation for the situations of figures 
(circles), as well as for the case |Q| > 1/2 (diamonds). The average radius of the drops is identified by the first 
zero of correlation functions. 



While the simulations confirm the scenario described in Sec. IV, we cannot extract the power law exponent predicted 
by the non conserved ordered parameter dynamics at short times, as well as those of the conserved order parameter 
dynamics for long times. Indeed finite size effects prevents us from reaching the conserved regime. 
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VI. CONCLUSIONS 



We have studied a model which reproduces many of the features which render appeaHng the scalar Phase field 
model and analysed its equilibrium and off-equilibrium properties. We transformed the original scalar model into a 
model with global couplings, more amenable to analytic investigations. 

In such a description the temperature shift from the coexistence temperature plays the role of an annealed field, 
which changes during the process and settles to a value determined by energetic considerations. Its dynamics is slow 
compared with that of the order parameter field and the latter becomes eventually slaved by the first. 

The long time state can be either a pure state with vanishing correlations, in the limit Tf = 0, or a mixed state 
with large spatial fluctuations in the order parameter. The type of equilibrium reached depends on the initial value 
of the spatial average Q of the field U. Indeed, in the case —1/2<Q< 1/2, the system shows a tendency towards 
separation into two phases in proportions given by the rule m = 2Q and one observes drops of the minority phase in 
a sea formed by the majority phase. At the same time the thermal field u vanishes, indicating that T = Tm in the 
whole volume. The number of the drops decreases with time as to minimize the free energy of the system, but for long 
times the total amount of solid remains fixed because the heat released by a growing solid drop can only be adsorbed 
by a shrinking solid drop. As a result the solid order parameter (j) becomes nearly conserved being mediated by the 
conserved heat field. At this stage the dynamics of the crystalline order parameter becomes a genuine conserved order 
parameter dynamics. One can in fact observe multiscaling, if the volume of the systems is large enough. The existence 
of inhomogeneous structures is mirrored in the presence of the peak in the structure factor at finite wavelength and 
the phenomenon is similar to the Ostwald ripening. The undercooling initially present is not sufficient to promote the 
transformation of all the liquid into the solid state, and some drops of either phase remain trapped into the other. 

To summarize the results we have established the following rules governing the evolution: 

1) The field U is constant in time. 

2) The order parameter m{t) as t — > oo tends to the asymptotic value m = 2Q, if Q falls in the range [—1/2, 1/2]. 
This fact in turn implies that the spatial average value of u over the system vanishes as t — > , i.e. the system reaches 
two phase coexistence asymptotically. 

3) In the above range of Q, the correlation function is large and centered at finite values of the wavevector k. 

4) If IQI exceeds the threshold value 1/2 the system evolves towards a spatially uniform state with m ^ —1 and 
vanishing correlations and m is no longer equal to 2Q. In this case u reaches an equilibrium value which is non zero 
and the system is out of two phase coexistence. 

5) If IQI > 1/2, i.e., larger values of the undercooling cause the melt to crystallize completely and thus correlations 
are asymptotically suppressed as the system reaches an homogeneous state. 

The above features are interesting because mimic the behavior of the more realistic Phase field model with local 
constraint. 

We, finally, remark that perhaps, the most serious flaw of the model, is that it suffers from the same problem as the 
spherical model. In contrast with the local constrained Phase field model, which displays a region of metastability of 
the liquid phase in the ^-T plane , between the coexistence line and the spinodal line [p6| , the globally constrained 
model is always unstable inside the two phase coexistence line. As a consequence, no nucleation barrier needs to 
be overcome in the transformation from liquid to solid. In the initial state long-wavelength fluctuations grow and 
the system becomes unstable. The nucleation barrier being proportional to the surface tension associated with the 
creation of a kink in the scalar model, whereas in the model with global couplings the energy gaps between the ordered 
phase and the instanton solutions, i.e. the uniform solutions of the equation vanishes in the infinite volume limit. 
Thus the mechanism described is non-Arrhenius like. 

This is also reflected in the absence of true phase separation, since the width of a domain wall diverges in the limit 
of a vanishing pinning field, h, as [ p3| . 
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FIG. 1. Schematic behaviour of dm/dt as a function of m, from eq. (p6|). 

FIG. 2. Typical behavious of the order parameter m as function of time for the cases: (a) Q > 1/2; (b) < Q < 1/2; (c) 
— 1/2 < Q < 0; (d) Q < —1/2. In all runs shown here we started with m(t = 0) = 0.54, while the other parameter were: 
(a) u{t = 0) = 0.4, Q = 0.67, m(oo) 2Q; (b) u{t = 0) = -0.2, Q = 0.07, m(oo) = 2Q; (c) u(t = 0) = -0.4, Q = -0.13, 
m(oo) = 2Q; (d) u(t = 0) = -1.0, Q = -0.73, m(oo) / 2Q. 

FIG. 3. The order parameter m as a function of time for Q — —0.3077 and different values of A, from rigth to left A — 128, 
170, 200, 300, 500. The plateau increases as A decreases. At the initial time the parameters are m = 0.385 and u — —0.5. The 
dashed line denotes the value 2Q. 

FIG. 4. The cross-over time tc as function of A for the curves of Fig. ^. The dashed line is the scaling 1/A. The cross-over 
time is defined as mint m{t) — 2Q. 
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FIG. 5. Snapshot of the 4> field for the global constraint case and —1/2 < Q < 0. On the axis we report the lattice index. 



FIG. 6. Snapshot of the 4> field for the local constraint case and — 1/2 < Q < 0. On the axis we report the lattice index. 

FIG. 7. Circularly averaged 4>4> correlation function as function of the lattice index for the global constraint case obtained 
numerically: circles represent \Q\ < 1/2; diamonds refer to \Q\ > 1/2. 

FIG. 8. Circularly averaged ipip eorrelation function as function of the lattice index for the local constraint case obtained 
numerically: circles refer to \Q\ < 1/2 whereas diamonds \Q\ > 1/2. 
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